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ABSTRACT 

The goal of the Flight Dynamics Facility (FDF) attitude support is twofold: to deter- 
mine spacecraft attitude and to explain deviations from nominal attitude behavior. 
Attitude determination often requires resolving contradictions in the sensor observa- 
tions. This may be accomplished by applying calibration corrections or by revising 
the observation models. After accounting for all known sources of error, solution 
accuracy should be limited only by observation and propagation noise. 

The second half of the goal is to explain why the attitude may not be as originally 
intended. Reasons for such deviations include sensor or actuator misalignments and 
control system performance. In these cases, the ability to explain the behavior 
should, in principle, be limited only by knowledge of the sensor and actuator data 
and external torques. 

This paper documents some results obtained to date in support of the Cosmic Back- 
ground Explorer (COBE). Advantages and shortcomings of the integrated attitude 
determination/sensor calibration software are discussed. Some preliminary attitude 
solutions using data from the Diffuse Infrared Background Experiment (DIRBE) 
instrument are presented and compared to solutions using Sun and Earth sensors. A 
dynamical model is constructed to illustrate the relative importance of the various 
sensor imperfections. This model also shows the connection between the high- and 
low-frequency attitude oscillations. 


*This work was supported by the National Aeronautics and Space Administration (NASA) /Goddard 
Space Flight Center (GSFC), Greenbelt, Maryland, Contract NAS 5-31500. 
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1. INTRODUCTION 


This paper compiles an assortment of results and comments regarding attitude determination and sensor 
calibration for the Cosmic Background Explorer (COBE). The central problem is to disentangle the true 
spacecraft motion from the apparent motion reported by imperfect sensors. When these sensors also 
control the spacecraft, their biases not only produce fictions in the measurements but induce wobbles in 
the true motion. 

First, a brief description of the COBE Attitude Control System (ACS) is given. The Diffuse Infrared 
Background Experiment (DIRBE) instrument is also described, in the context of using star sightings to 
increase attitude determination precision. Section 2 discusses the advantages and shortcomings of the 
integrated attitude determination/sensor calibration software, with emphasis on real-world difficulties. 
This is followed by a comparison of preliminary attitude solutions using DIRBE data with results using Sun 
and Earth sensor data. The final section is meant to complement the data reduction techniques of Sec- 
tion 2. Rather than solving for the sensor biases that best match the predicted observations to the real 
data, the dynamical equations themselves are solved. This allows the separate study of the effects of each 
bias or misalignment. 

1.1 DESCRIPTION OF THE CONTROL SYSTEM 

COBE is a spinning, three-axis stabilized spacecraft. Two counter-spinning momentum wheels control the 
spin rate and leave the spacecraft with zero net angular momentum (except for the 1 -rotation per orbit 
(rpo) pitch rate about the Sun line) . Sensors, actuators, and much of the control electronics exist in 
triplicate on the three control axes (A, B, and C) . These axes are symmetrically located in the plane 
normal to the spin axis. The sensors include a gyroscope, an Earth scanner assembly (ESA), a digital Sun 
sensor (DSS), and a three-axis magnetometer (TAM) on each control axis. The actuators are reaction 
wheel assemblies (RWAs) and magnetic torque rods (transverse and X-axis) for momentum unloading. 
The ACS drives the three reaction wheels independently on each axis in proportion to rate, roll, and pitch 
error signals (Reference 1). 

The DSS most directly facing the Sun measures the Sun elevation and azimuth. A common electronics 
unit produces the sine and cosine of the Sun azimuth angle. It also produces the sine and cosine of this 
angle plus or minus 120 degrees (deg) for use by the other axes. 

The rate error signal is the gyro output minus an orbit rate-stripping term, normally equal to 1 rpo multi- 
plied by the cosine of the Sun azimuth. Choosing a rate-stripping parameter different from 1 rpo gener- 
ates an error signal that can be balanced only by a nonzero pitch signal. For this reason, the rate-stripping 
parameter also serves as the pitchback parameter. 

Roll angle is minus the Sun elevation angle. The commanded roll offset is subtracted, and the result is 
multiplied by the sine of the appropriate Sun azimuth to yield the roll error signal. 

The pitch signal is the difference between the instantaneous ESA split-to-index angle and a modulated 
reference value. The reference is the split-to-index angle sampled when the Sun azimuth is 90 deg (actu- 
ally 81 deg to avoid Sun interference) . At 90 deg, the ESA is directly reading the roll angle. This value is 
held for a full spin period. It is multiplied by the sine of the Sun azimuth and subtracted from the 
instantaneous split-to-index. With the roll contribution thus stripped out, the result is the pitch angle 
multiplied by the cosine of the Sun azimuth. 

The rate, roll, and pitch error signals are amplified, filtered, limited, and combined into a torque com- 
mand signal for the RWAs. The reaction wheel tachometer also provides feedback so that, in the absence 
of other signals, each wheel will spin down to a commanded speed. 

The DSS has a precision of only 0.5 deg, which is an inherent limitation on the ability of the ACS to 
establish the roll angle and the correct azimuthal phase for modulating the sensor output. Star sightings in 
the DIRBE instrument, described in the next section, offer data with greater precision. Although these 
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measurements are not available to the ACS for attitude control, they can be used a posteriori for im- 
proved ground-based attitude determination. 

1-2 DJRBE STAR IDENTIFICATION 

DIRBE is a COBE science instrument whose mission is to perform a full-sky survey of diffuse infrared 
radiation in the wavelength range of 1 to 300 microns. As a bonus, DIRBE also provides star sightings to 
the Flight Dynamics Facility (FDF) that enhance the FDF’s ability to determine attitude. DIRBE has a 
square field of view, 0.7 by 0.7 deg. The DIRBE boresight (at the center of the field of view) is oriented 
30 deg from the spacecraft spin axis. The field of view sweeps out a spiral pattern on the celestial sphere 
as COBE spins. Two edges of the square field of view are parallel to the scan direction. 


As DIRBE scans the sky, pointlike sources, such as stars and planets, pass through the field of view. 
These sources lead to a sharp increase in the measured infrared intensity. Examples of these “spikes” in 
intensity are shown in Figure 1. The spike profile can be processed (Reference 2) to determine when the 
pointlike source passed through the center of the field of view. This “time of passage” (TOP) is then used 
to predict the boresight direction by interpolating approximate attitudes computed from the attitude deter- 
mination subsystems. 

The predicted boresight direction at the TOP is used for star identification. The observed star is com- 
pared to the SKYMAP 2.2 micron (K-band) wavelength reference catalog (Reference 3). For a positive 
star identification, a reference star must be located within a user-specified angular radius of the predicted 
DIRBE boresight. Approximately 2000 well separated, sufficiently bright stars are currently in the refer- 
ence catalog, and about 1.3 valid DIRBE star identifications occur per minute. 

Individual star sighting accuracies are expected to be ±0.35 deg in COBE body elevation and ±0. 1 deg in 
COBE body azimuth. These star sightings have been used to determine the attitude of COBE independ- 
ently of Sun and Earth data. In turn, these attitude solutions can be used for sensor calibration. 


DIRBE D*!A VS t I ME 



Figure 1. Converted DIRBE Data 
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2. DETERMINING THE ATTITUDE 


Two of the tools used to determine the attitude of COBE are the Coarse Attitude Determination Subsys- 
tem (CADS) and the Fine Attitude Determination Subsystem (FADS). 

CADS computes single-frame attitudes from “simultaneous” sensor measurements that have been interpo- 
lated to a common time. As the content of ESA and DSS data is partly redundant, the user can choose 
how heavily to weight each. This choice results in a range of solutions reflecting the different sensitivities 
and misalignments of the various sensors. 

FADS estimates these misalignments, the gyro biases, and an epoch attitude by performing a least-squares 
fit to the data over a given timespan. The bias-corrected gyro data are used to propagate the epoch 
attitude over the timespan. 

2.1 INTEGRATED ESTIMATION 

2.1.1 BACKGROUND 

COBE is the first spacecraft supported by the FDF to estimate attitude, gyro, and sensor calibrations in a 
single ground support subsystem, FADS (Reference 4). This integrated capability is the next logical step 
from the Earth Radiation Budget Satellite (ERBS) FADS (Reference 5), which combined attitude and 
gyro calibration, and follows in the spirit of the Multisatellite Attitude Determination/Optical Aspect Bias 
Determination (MSAD/OABIAS) System (Reference 6), which is the standard estimator for spinning 

spacecraft. 

The COBE FADS is a batch-weighted least squares estimator that can use any combination of the 23 ob- 
servation types to solve for almost any combination of the 49 possible state parameters. The three attitude 
parameters must be solved for as a group. This approach is attractive because it obviates the need for 
separate calibration utilities and the need to go back and forth between them to obtain a complete set of 
self-consistent calibrations. Instead, all parameters are solved for in the same subsystem and, in principle, 
at the same time. Although the first benefit is definitely realized, the second benefit, simultaneous solu- 
tion, is problematic for several reasons. 

First, solving for many parameters at once often causes the solution to diverge or to converge to an answer 
worse than the a priori solution. Only when the parameters being solved for are very close to their correct 
values is there a benefit to solving for many parameters at once. Then, the solution usually improves. 
This problem may be solved, in whole or part, by changing the method of solving the batch least-squares 
problem. Marquardt’s algorithm and others offer alternatives to the standard Newton-Raphson met od 
used here. 

Second, a byproduct of this great flexibility is confusion. With so many parameters, it becomes hard to 
decide what to include in the state vector. At present, choosing parameters depends on the analyst s 
ability to recognize patterns in the residuals and to identify the parameter causing the pattern. Because 
many parameters have similar effects and may combine to produce unrecognizable residual patterns, it 
may take longer than the 2 months usually allotted for launch support to decide exactly what is wrong. 

Third, whether it is a simple mistake in the specifications or a more fundamental problem, gyro calibration 
takes a very longtime to converge. As COBE convergence criteria are maximum values arbitrarily set for 
the change to the state, a more meaningful comment might be that the gyro calibrations approach their 
final values slowly and monotonically. They never overshoot their true values, as do the other parameters. 
Moreover, they approach slowly, typically by halves. Because FADS is already slow, the need for so many 
iterations increases the central processing unit (CPU) time demands and restricts calibration wor to 
off-hours when the mainframe is not being heavily used. 

Thus, although it is possible to solve for all 49 parameters at once, that has not been a real benefit of the 
integrated calibration approach. The main benefits have been in centralizing the solution process and in 


82 



providing a better, more complete view of the calibration problem, as reflected in part by the following 


2.1.2 STATISTICS 

Statistics are useful for assessing solution accuracy and selecting parameters for the state vector. The 
following statistics are being tried for COBE. Knowing the solution accuracy can help in recognizing when 
a solution is diverging and in choosing between competing solutions. To evaluate the quality of a batch 
solution, the FADS provides the weighted root mean square (RMS) residual (a,). This is a single number 
whose magnitude is the size of a typical residual (A yi ) at time t s weighted by the weight matrix (W). The 
set of calibrations that gives the smaller RMS for the same data is probably the better set. To first order 
the RMS residual is independent of the number of data points and the choice of observation weights. 


<% = 




AyT W Ay s l/N trace (W) 


( 1 ) 


To know how accurate each parameter of the solution really is, FADS provides a residuals-based estimate 
°. * ® covariance matrix. Rather than assume that the observation residuals are zero mean and white, 
with the variance assumed for weighting purposes, the covariance estimate is computed from its definition: 


cov (x) = E [Ax Ax t ] - E [Ax] E [Ax] T (2) 


The two terms are computed from the expressions for the change to the state. Here, F, represents the 
partial derivative of the observation types with respect to the state variables at time i, based on the same 
currently estimated state used to predict the observations.* 



One possible explanation for the divergence of solutions and the slow convergence of the gyro calibration 
is that the observabilities of the state parameters are too different. For a state vector of two elements, this 
means that the level curves of the loss function being minimized are elliptical rather than circular. If the 
minimization algorithm moves perpendicular to these level curves, the less observable parameter can be 
slow to converge (Figure 2). 
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Figure 2. Convergence With Very Different Observabilities 


The diagonal elements of the normal matrix (2 F T W F) provide estimates of the relative observability of 
the different parameters in the state vector. Due to differences in the units used for computation and 
display, however, the standard deviations in the iteration summaries may not reflect the true conditioning 
of the normal matrix. A more direct indication of the unbalanced observabilities and singularity is the 
condition number (K), which is the product of the row norm of the normal matrix and that of its inverse 

(Reference 7): 


K(A) = || A | | A' 1 || 


When K is large, the parameters solved for have very different observabilities, and the matrix is nearly 
singular. As a further check, the maximum and minimum eigenvalues of the normal matrix are computed 
along with their eigenvectors. The associated eigenvectors indicate which parameters are the most or least 
observable. This tool has not proven very useful because condition numbers routinely reach 1 million 
without the solution diverging, and convergence can be slow without there being a very large condition 
number. 

Knowing what is most observable and least observable, however, does not dictate the parameters to in- 
clude in the state vector. The most observable parameter may be completely accurate, and the error may 
be caused by the least observable one. An enhancement, which has not yet been implemented, takes the 
inner product of the residual histories over the batch timespan with the predicted residual histories due to 
each parameter error. The computation entails normalizing the derivative of the residual vector (Ay) at 
time i with respect to each parameter ( x s ) : 



( 5 ) 


where 


k 


dAy t 

dxj 
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and accumulating this dot product to indicate which parameters best fit the residuals. 


N 

^ A Vi 


( 6 ) 


i = 1 


2.1.3 COBE CALIBRATION HISTORY 

Before launch, calibration looked straightforward. FADS was working. Signatures of the individual cali- 
brauon errors had been identified, and a “foolproof” procedure had been set down whereby anyone 
could determine the calibrations. The only thing that would prevent calibration from being completed the 
day of launch was that DIRBE measurements would not yet be available. 

Important results were obtained in the first few days. Earth sensor-A appeared to have a twist misalign- 
ment (spht-to-index bias) of 0.5 deg, and the X-gyro scale factor was 0.1 percent high. After that, 
however, progress was slow. Finding scale factor corrections for the transverse axis gyros proved treacher- 
ous, and when transverse gyro biases were solved for, the results were not repeatable. Localized sources 
o error, such as horizon radiance, also turned out to be more significant than expected. 

One of the early “casualties” in this struggle was the original calibration plan. Reasonably, it had seemed 
sensor alignments should be solved first, over short timespans, before gyro errors could become signifi- 
cant. Gyros would then be calibrated. In practice, however, the X-gyro scale factor error could not be 

ignored even over a single 75-second spin period, and the procedure was changed to make it the first 
parameter solved. 

With the Earth sensor and X-gyro scale factor corrections, it was possible to increase the length of the 
batches to 20 minutes without the observation residuals exceeding 0.5 deg. Prelaunch information had 
also suggested that batches longer than this would have propagation errors of more than 0.1 deg due to 
gyro noise. As a result, 20 minutes became the standard batch length. The next significant problem was 
the growing Earth sensor split-to-index residuals. These residuals could be reduced by assuming large 
transverse scale factor errors. This solution was considered valid for several weeks until it was noted that 
attitude histories computed in this way had pitch discontinuities at the batch boundaries. The large scale 
factors had to be abandoned. 


After a few weeks, it was realized that without the large scale factors the roll drifted off, and solving for 
sensor alignments along with the X-gyro scale factor made the residuals and discontinuity much smaller 
When a full orbit was processed, the roll came back over the second half of the batch. This identified the 
problem as a spin phase error similar to that caused by incorrect telemetry time tags or single-step attitude 
propagation (Reference 8). It turned out that the Sun sensors that measure yaw were slightly misaligned 
with respect to the gyros. By correcting the Sun sensor azimuth alignments by approximately 0 1 degree 
the problem was resolved. 


Recently, attention has returned to the question of transverse gyro scale factor errors, because single- 
frame and FADS pitch solutions diverged slightly over long times. When 0.002 scale factor corrections 
were made, the Earth sensor residuals were halved, and the gyro-propagated attitude repeated from orbit 
to orbit. These results corroborate earlier ones obtained without examination of pitch repeatability, the 
full importance of which was not appreciated at the time. 


One of the first calibration objectives following launch was to determine transverse gyro biases. Even 
today, these remain a puzzle. The problem is that Earth sensor twist misalignments and transverse gyro 
biases have the same effect on the split-to-index residuals. If zero biases are assumed, Earth sensor-A 
twist appears to be 0.5 deg. If the twist is fixed at its prelaunch value of 0.3 deg, the biases on gyro- A and 
-C are on the order of 30 deg per hour. It has not been determined what the real answer is. 
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The small coning motion implied by such a twist or bias is not readily observable from the Sun sensor 
elevation measurements because of their coarse 0.5 deg quantization and the attitude being controlled to 
maintain constant Sun elevation. DIRBE elevation measurements, with their 0.7 deg field of view and 
approximate 1-minute spacing, are similarly crude. Nor does the control system provide any help. As 
both the Earth sensors and the gyros are in the control loop, various mixtures of the two errors also cause 
coning. This may be a situation in which if it is so hard to decide, it probably does not matter. 

Finally, although observation residuals are now about as small for a full orbit as they are for one spin 
cycle, old data must be reexamined to see if the calibrations that work well now would have been as 
successful then. Horizon radiance effects, which peak at the solstices and vanish at the equinoxes, might 
have contributed to those large early residuals. Other localized errors that appear at the orbit rate and its 

harmonics are also under study. 

The COBE calibration effort has not been quick or smooth; each advance has been hard won. Moreover 
these solutions could probably not have been foreseen before launch, even with twice the preparatio 
time. The combined spinning and pitching made COBE calibration more confusing, but each mission is 
certain to have its own challenges, and so calibration may always be a puzzle. In solving the puzzle, ideas 
are like pieces, the more the better, but each must be evaluated and recorded to avoid confusion and 
needless duplication of effort. The integrated estimation software was able to reduce some of this confu 
sion by providing a broad view of all the observations and calibrations at once, but it was not a panacea. 

2.2 DIRBE ATTITUDE RESULTS 

Preliminary results with DIRBE indicate good agreement with Sun and Earth attitude solutions. However, 
DIRBE misalignments have not yet been computed, and the solutions seem sensitive to the^ misa- 
ignments. A comparison was made between FADS runs using data from January 29, ■ ( 

nominal case) assumed the DIRBE boresight to be aligned with the nominal body azimuth (240 deg in the 
body). The other (the misaligned case) assumed a -0.3 deg misalignment of the boresight with respec o 
the nominal azimuth. The computed roll in the nominal case, when averaged over a ^' nC ^ le “ e ™'" 
nate the effects of short-term oscillations, experienced a long-term drift of 0.3 deg over 25 minutes ( g 
ure 3). This agreed with FADS solutions, using Sun and Earth data, to better than 0.06 deg m the roll, 
pitch, and yaw. In the misaligned case, however, roll experienced a long-term drift of only 0.1 deg in 
25 minutes (Figure 4). Sun elevation data drift by roughly 0.1 deg over the 25-minute time interval. 
Based on Sun elevation information alone, the misaligned DIRBE solution is more believable. However, 
the half-degree digitization of the DSS makes such a comparison suspect. Adding to the uncertainty is the 
possibility that the Earth scanner information is corrupted by unmodeled horizon radiance effects. 


3. DYNAMICAL MODEL 


This section approaches the attitude problem from a different direction. The intent is to understand what 
is driving both the high-frequency motion (time scale = spin period) and the low-frequency motion (time 
scale = orbital period), which have been observed since the failure of the B-axis control gyro. The starting 
point is the system of Euler equations and the attitude control laws. Sensor misalignments, biases, and 
failures can be added one at a time to see how each contributes to the attitude motion. 


During one spin period, the COBE body X c -axis traces out two unequal cones; the pitch and roll angles 
oscillate with dominant frequency components equal to the spin rate and twice the spin rate. e amp i 
tudes are a few tenths of one degree and are dependent on the spin rate. 


During one orbit period, the 
pending on the time of day. 
nents are (o and 2 o ) , where 


spin-averaged pitch angle varies by a few degrees, with the amplitude de- 
Two unequal peaks appear each orbit, so the dominant frequency compo- 
o) is the orbital frequency. 
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Figure 3. Roll Angle Solution Using Figure 4. Roll Angle Solution Using -0.3 Deg 

Nominal DIRBE Alignment Azimuthal DIRBE Misalignment 

3.1 MODEL EQUATIONS 

The derivation of the equations of motion and control laws is outlined in Appendix F of Reference 1. 
The following is a brief summary of those equations and definitions. (Reference 9 also derives similar 
equations and applies them to a stability analysis of the nominal ACS configuration.) 

A convenient reference system is the despun frame, which is the COBE body frame rotated to zero deg 
yaw angle. The Euler equations in this frame can be written 


L x = Q cos 9 h r - 0 h p + T x 


10 = - Q cos 0 L x - I Q 2 sin 6 cos 9 - Q sin 0 h p - h r + T r (?) 


I (<P ~ V) cos 0 = 2 I Q 0 sin 9 + 9 L x + £2 sin 9 h r - h p + T p 


where 


V'* 0» <t> 

= yaw, roll, pitch 

V 

= ideal pitch anomaly 

Q 

= (p - V + a> 

(!) 

= orbital rate 


I = transverse moment of inertia = 2700 ft-lb-sec 2 


NAL PAGE !S 
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L x = X-axis angular momentum 

h r , h p = net angular momentum of the three reaction wheels, measured in the body 

frame and transformed to the despun frame 

T x , T r , T p = external torques (e.g., gravity gradient and torquer bars) 

The ACS control torque commands to the three wheels are projected on the body Y c - and Z c -axes and 
then transformed to the despun frame: 

K t h r + 1.5 (K p + K r F 20 ) {0 - 0 O ) + 1.5 K g 0 - h p 

1 °) 

K, h p + 1.5 K p 0 + 1.5 Kg (Q cos 0 - Q C md) + V’ h r 

where 0 O is the roll offset and Q cm d is the orbit rate-stripping parameter. The feedback loop gains are 

K p = 1.07 ft-lb/rad 
K r = 6.59 ft-lb/rad 
Kg = 189.6 ft -lb/ (rad/sec) 

K t = 0.0168 sec' 1 



The yaw rate is related to the X-angular momentum through the following equation: 

ip = ip 0 + L x /I x + sin 0 


( 9 ) 


where ip 0 is the nominal yaw rate. 


The ACS roll loop has a 20-second filter to smooth out the 0.5-deg digitization of the DSSs. 
digitization is not included in the model, the 20-second filter has been retained. The filter 
converts 0 into Or where 


Although the 
operator F 20 


(20 sec) 0 f + 0f = 0 - 0o 


To a good approximation, the pitch anomaly is 




360° 


tan 2 (<p/2) sin (2cot) 


( 10 ) 


(ID 


where (p is the difference between the Sun’s declination and that of orbit normal. 
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Linearizing the dynamical system in the roll, pitch, L x , and control torques, and taking the Laplace 
transform, yields 



where s is the complex Laplace transform variable, A <p is the transform of <j> — tj , and 


A = | (K p + s Kg) 

B = A ♦ | K,/(l ♦ 20s) 

E = K, «„■ ( *1 

2 p \ 2 / 4o> 2 + 

P = | Kg (£2 - Q cmd )/ S 


(13) 


The initial conditions have all been taken to be zero. P is the pitchback command, and E is the pitch 
anomaly. H r , H p are errors caused by the failed gyro and ESA twist misalignment discussed in the next 
subsections. T x is assumed to be zero. 

The characteristic equation is obtained by setting the determinant of the matrix in Equation (12) to zero. 
The determinant can be written 


det = I 2 (s 2 + w 2 ) D 


(14) 


D = 


1 

0.05 + s 


n (* - Sk) 


ksl 


(15) 


89 

r \ 



With V > 0 = -0.815 revolutions per minute (rpm) and two working gyros, the roots are 


si = -0.009715 

s 2 = -0.02740 + i 0.01604 ( 16 ) 

s 3 = -0.07978 + i 0.08095 


and the complex conjugates, all having units of sec' 1 . These are the time constants for the decay of 
transients. They all have negative real parts, so the system is linearly stable. 

One can solve the algebraic equations for and then obtain A0(t) by inverting the Laplace transform. 
The inverse is obtained by inspection from the singularities in A<p, discarding the transient solutions. The 
residues of the simple poles give the amplitudes. Double poles lead to oscillations plus secular terms; that 
is, terms that grow linearly with time. These occur when an external roll torque oscillates at the orbital 
rate, always pointing the same direction in inertial space. Such secular terms (for example, the gravity 
gradient torque) must be canceled by the secular terms from the momentum management assembly 
(MMA) torque. On average, the spacecraft is not accumulating angular momentum. These terms are 
dropped in the following solutions. 

In this manner, the equations have been solved for the commanded pitchback (PB) , the pitch anomaly 
(PA), and the gravity gradient (GG) source terms; 


<p(t) = -6.883 (PB) 

-0.1833 sin (2<ut + 76.75 deg) (PA) 
-0.3721 sin (rut - 34.14 deg) (GG) 


(17) 


in units of degrees. These solutions assume that the yaw rate is -0.815 rpm, the Sun declination is 
-21 deg, the roll offset is -4 deg, and the gain in the gyro feedback loop is appropriate for two working 
gyros. The pitch anomaly has been added back in, so Equation (17) represents the actual pitch angle and 
not A <p . 

3.2 B-AXIS GYRO FAILURE 

With the failure of the B-gyro, the ACS torque command to the B-axis RWA became 


h B = hB nom) - K g [6 sin Vb + (0 - 9 + AQ ) cos 
AQ — 0) — £2 C md t lift 


(18) 


where h B nom) is the nominal torque command, V'b = - (t/> + 120 deg) is the Sun azimuth measured 
from the B-axis (note that azimuth and yaw angles are opposite in sign), Gcmd is the rate-stripping pa- 
rameter on the A- and C-axes, and £2 B is the rate-stripping parameter for the B-axis. Since the gyro 
failure, rate-stripping has been set to zero on all three axes, but it may be desirable to command a 
different value for the B-axis now that the symmetry has been broken. To allow for this, the nominal 
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£2 C md is subtracted, and £2 b is explicitly included. The sine and cosine factors in Equation (18) repre- 
sent the projection of roll and pitch rates on the B-axis. 

After projecting the new B-axis torque command on the body Y c - and Z c -axes, transforming to the 
despun frame, and taking the Laplace transform, the resulting source terms become 


H r = 



+ 



+ L [d cos (2y> + 


60°) + A (j> sin (2 ip + 60° 



(19) 


H p + P = 



+ L [0 sin (2t{> + 60°) - A0 cos (2^ + 60°)] 


} 


( 20 ) 


The first term on the right-hand side of Equations (19) and (20) can be moved to the left-hand side of 
Equation (12). This effectively reduces the rate gain from 3 Kg/2 to K g in the matrix elements A and B. 
The second term in both equations drives the system at twice the spin rate. 

The third term in Equation (20) determines the pitchback angle when it is balanced by the pitch feedback 
term from Equation (8). The nominal pitchback term, P, is included in Equation (20) to make explicit 
the net pitchback. Writing 2Q cmd as £2 a + Qc, the pitchback angle for the ACS with N„ working 
control gyros is 


Kg ( Ng o) Q a + Q b + Q c 

K p ^ 3 3 


( 21 ) 


With the rate-stripping parameters set to zero and N g = 2, the pitchback is -6.883 deg. 

In the final terms of Equations (19) and (20), the L is the Laplace transform operator. These terms are 
dropped in first order. They are added back in iteratively, using the rates obtained in the first-order 
solution. This produces terms with zero frequency (constant 0.1-deg roll and 0.75-deg pitch offsets) and 
small oscillatory terms with a frequency four times the spin rate (amplitude 10~ 5 deg). 

In solving for roll and pitch as functions of time, it is found that the failed gyro forces the spacecraft to 
move on an almost circular cone, twice per spin: 


6 = 0.100 - 0.103 sin (2ip + 14.46°) 


< P — -6.883 + 0.745 — 0.102 sin (2 y> - 76.60°) 


(23) 


where ip - ip a t and H> 0 = -0.815 rpm. Angles are given in degrees. The -6.883 deg term is the normal 
pitchback, and the other two constants are the second-order corrections. The cone is larger and less 
circular for lower spin rates. 
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The part of the high-frequency motion that is driven by the failed gyro can be prevented by choosing 
rate-stripping parameters that make the oscillatory source terms vanish; that is, set Aft = 0 . This can be 
done for any desired pitchback angle, <l>o. The required parameters are 

fi A = Qc = ^cmd 

ft B = Gcmd - a> ( 24 ) 

ftcmd = 0 > + <Po (Kp/Kg) 


The choice of ft A and ftc is not affected by the failed gyro, and the value of ft B is reduced by <u. 

3.3 ESA TWIST M ISALIGNMENT 

The sensor calibration feature in FADS has established the ESA twist misalignment angles to be T A = 
0.54 deg and T c = -0.15 deg. These twists alter the index angle in measurements of the spht-to-index. 
(This FADS solution assumes that all the error is due to ESA twist and none to transverse gyro drift bias.) 

Analysis similar to that of the previous section yields the attitude solutions (in deg) : 

6 = -0.0182 - 0.0394 sin ( ip - 85.38°) + 0.0063 sin (2ip + 32.54°) 


(25) 


<j> = +0.0355 sin (ip + 0.72°) + 0.0062 sin (lip - 58.53°) 


(26) 


The error in the instantaneous split-to-index drives the system at the spin rate; the error in the sample- 
and-hold split-to-index drives the system at twice the spin rate and generates the offset. 

The predicted high-frequency attitude motion is the sum of these expressions and the failed gyro solutions 
given in Equations (22) and (23). This motion is plotted in Figure 5. The figure shows two spin periods, 
starting from ip = 0 at t = 0. Figure 6 shows two spin periods of spacecraft data as determined by 
FADS. In that figure, the yaw angle is zero at time 19.31.11. The phase and overall shape of the FADS 
solution are well reproduced by the model equations driven only by gyro and ESA errors. The average 
values disagree because neither the commanded roll offset nor the magnetic torques are included in 
Figure 5. The predicted amplitude falls below the actual value by approximately 0.05 deg. These calcu- 
lations have been repeated for % = -0.225. -0.4, and -0.6 rpm. The qualitative agreement with FADS 
is good at all spin rates, but the amplitude is consistently low. This is possibly caused by an unmodeled 
momentum wheel misalignment, which also shows up as a larger-than-predicted, spin-rate-dependent off- 
set of the reaction wheel speeds. Gyro drift biases of roughly 30 deg per hour also could possibly make up 
the difference. 

3.4 MAGNETIC TORQUES 

Under ideal conditions, the RWA holds the spacecraft at the desired attitude with the wheels absorbing all 
the accumulated angular momentum. The MMA takes the tachometer signals, filters out any constant 
offset, and energizes the torquer bars to dump this excess angular momentum. However, imperfect sen- 
sors cause the ACS to spin the wheels up and down, even when the attitude is nominal. The MMA 
receives the oscillating wheel speed signals and attempts to dump excess momentum when there is none. 

The model assumes a dipole Earth field, oriented so that COBE’s orbit passes over the dipole axis at the 
northernmost point at t = 0. The field components in the body frame depend primarily on the yaw angle. 
For simplicity, they are calculated with pitch and roll assumed as zero. The reaction wheel speeds are 
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Figure 5. Roll and Pitch Angles Driven by ESA Misalignment and Failed B-Gyro 
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Figure 6. FADS Roll and Pitch Solutions 
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found by integrating the ACS control laws, using the high-frequency attitude solutions already obtained. 
In agreement with spacecraft data, the A-, B-, AX-, and BX-torquers are driven to their limits by the large 
oscillations of the A- and B-wheels. 

Transforming the A-, B-, C-dipoles to the despun frame multiplies them by sine and cosine of the yaw 
angle. This beats with the reaction wheel frequencies (predominantly the spin rate) . As a result, the pitch 
and roll components of the MMA dipole have both an offset and an oscillation at twice the spin rate. The 
MMA flips the sign in the Southern hemisphere. Multiplying by the X c -component of the Earth’s mag- 
netic field and spin-averaging yields torques that vary as the absolute value of cos (a>t). 

The final step is to put the MMA torques into Equation (12) and solve for the attitude. The spin- 
averaged torques are crudely approximated to be 

T r = -0.9 x 10~ 3 |cos art| (ft -lb) (27) 


T p = -0.2 x 10' 3 - 1.1 x 10~ 3 | cos <nt| (ft-lb) (28) 


Solving for Atf>, one finds terms with single and double poles at s =0, and single poles at s = ±ift> and s - 
+i2n to (n running from one to infinity). The amplitudes drop off rapidly with n; only the n = 1 terms are 
kept. They drive the system at twice the orbital frequency. 

Figure 7 shows two orbits of the pitch motion driven by pitch anomaly, gravity gradient, and MMA 
torques. The spacecraft data are shown in Figure 8. The agreement is not good, although the model does 
give the large amplitude and two unequal peaks. The shape and amplitude are fairly sensitive to the 
parameter used in Equation (27), which in turn is sensitive to the details of the model (for example, 
inclusion of ESA misalignments and second-order gyro effects) . 

The solution could be improved by using a more realistic magnetic field model and by taking proper 
account of the attitude when calculating the field components in the body frame. Using the true attitude 
rather than the nominal values would also improve the gravity gradient results. This could be important 
because the gravity gradient natural period is approximately 93 minutes, comparable to the orbital period. 

With the absence of any other large perturbations, it is concluded that the MMA is the immediate cause 
of the low-frequency pitch motion. In turn, the MMA is responding primarily to the inconsistent 
pitchback signals due to the failed B-gyro. It is suggested that commanding the same pitchback angle on 
all three axes, as in Equation (24), could significantly improve spacecraft performance. 
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Figure 7. Pitch Motion Driven by Pitch Anomaly, Gravity Gradient, and MMA Torques 
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Figure 8. Two-Orbit CADS Pitch Solution 
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